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ABSTRACT 

The Coma cluster of galaxies hosts the brightest radio halo known and has therefore been the target 
of numerous searches for associated inverse Compton (IC) emission, particularly at hard X-ray energies 
where the IC signal must eventually dominate over thermal emission. The most recent search with the 
Suzaku Hard X-ray Detector (HXD) failed to confirm previous IC detections with RXTE and Beppo- 
SAX, instead setting an upper limit 2.5 times below their nonthermal flux. However, this discrepancy 
can be resolved if the IC emission is very extended, beyond the scale of the cluster radio halo. Using 
reconstructed sky images from the 58-month Swift BAT all sky survey, the feasibility of such a 
solution is investigated. Building on Renaud et al., we test and implement a method for extracting 
the fluxes of extended sources, assuming specified spatial distributions. BAT spectra are jointly fit 
with an XMM-Newton EPIC-pn spectrum derived from mosaic observations. We find no evidence for 
large-scale IC emission at the level expected from the previously detected nonthermal fiuxes. For all 
nonthermal spatial distributions considered, which span the gamut of physically reasonable IC models, 
we determine upper limits for which the largest (most conservative) limit is < 4.2 x 10~^^ erg s~^ 
cm~^ (20-80 keV), which corresponds to a lower limit on the magnetic field B > 0.2 /iC A nominal 
flux upper limit of < 2.7 x 10^^^ erg s^^ cm^^, with corresponding B > 0.25 /iC, is derived for the 
most probable IC distribution given the size of the radio halo and likely magnetic field radial profile. 
Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Coma) — intergalac- 

tic medium — magnetic fields — radiation mechanisms: non-thermal — X-rays: 

galaxies: clusters 
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1. INTRODUCTION 

The X-ray emission from clusters of galaxies is primar- 
ily thermal in origin and is produced by a diffuse popula- 
tion of intergalactic electrons in the ionized intracluster 
medium (ICM). These electrons coexist with a nonther- 
mal, relativistic electron population in at least some clus- 
ters - inferred from observations in the radio regime - 
which should also radiate at X-ray energies. While ther- 
mal emission clearly dominates in the kilo-electron volt 
(keV) energy range, it declines rapidly outside this range, 
allowing the detection of a nonthermal spectral signature 
as soft or hard excess emission. This possibility is espe- 
cially promising at hard (>10 keV) energies, where the 
exponential decline of the thermal bremsstrahlung con- 
tinuum is distinctly steeper than the expected nonther- 
mal spectrum. Measurements of nonthermal X-ray emis- 
sion are critical to the determination of the total amount 
of relativistic energy in the ICM, which is currently 
poorly constrained. While no more than ~10% of this 
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energy is tied up in nonthermal components, amounts at 
or near this level will af fect the dynamics and structure 
of the thermal gas (e.g-. IVazza et al.ll2009f ). Specifically, 
studies that attempt to infer the total masses of clus- 
ters from the hydrostatic state of the thermal gas will 
produce biased mass estimates if the pressure support of 
relativistic particles and fields is not accurately included. 
The mass functions built from these estimates can be 
used to constrain cosmological parameters; these studies 
are already und erwa y using observables derived in both 
the X-ray (e.g., M antz et al. 2008; Vikhhnin et al. 20091) 
and microwave (e.g., Va nderlinde et al.l l2010h through 
the Sunyaev-Zel'dovich effect) regimes. 

A measurement of the total energy in relativistic ICM 
components is possible when X-ray and radio nonthermal 
fluxes are combined. Diffuse, cluster- wide synchrotron 
radio emission, called radio halos or relics depending 
on their morphology, imply that both magnetic fields 
and relativistic electron populations are present on large 
scales. The total luminosity of a synchrotron-emitting 
electron is given by 



Lr = -cttct cb , 



(1) 



where cft is the Thomson cross-section, c is the speed 
of light, 7 is the Lorentz factor of the electron, and 
€b — B'^/8n is the energy density of the magnetic field. 
For a collection of relativistic electrons, the value oi Lr 
depends both on the number of electrons and on B and 
cannot independently determine either. However, these 
same electrons will up-scatter cosmic microwave back- 
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ground (CMB) photons through inverse Compton (IC) 
interactions, which have a luminosity Lx equivalent in 
form to equation ^ but with e b replaced by the energy 
density of the CMB. Since both luminosities are propor- 
tional to the number of electrons, their ratio gives the 
volume-averaged magnetic field, 



L 



X 



(2) 



CMB 



where a is the radiation constant and Tcmb is the tem- 
perature of the CMB. The IC r adiation should be ob- 
servable at hard X-ray energies (|Rephaelil[l977l) . Thus 
far, IC emission has only been detected at low sig nifi- 
cance ()Nevalainen et al.ll200l iMillion fc A llen"2009f) or, 
in o ne case at higher s ignificance, i n the Ophiuchus clus- 
ter (lEckert ct al. 2008': but see also Aiello ct al. 2009 and 
iFuiita etal. .2008) , although the diffuse radio emission 
in this cluster is res tricted to a smaller scale mini-halo 
(jMurgia et al.|[2010t ). The measurement of an IC flux 
from a synchrotron source directly leads to a simulta- 
neous determination of the average value of B and the 
relativistic ele ctron density (Harris & Romanishin 1974; 
iSarazinI 119881 ). Therefore searches for IC emission co- 
incident with a radio halo or relic are an excellent way 
to constrain the contribution of relativistic materials in 
clusters. 

The fi rst, and brigh test, radio halo was discov- 
ered by IWillsonI (jl970H in the Coma cluster, and 
its radio pr operties have perhap s been the best 
studied (e.g. iGiovannini et alj 119931: iDeiss et all 119971: 
iThierbach et al.l I2003D . Coma has been observed by 
all the major observatories with hard X-ray capa- 
bilities (Henriksen & Mushotzkv 1986"; 'Bazzan o et al.l 
[1990; Hughes et al. 1993; Rcphaeli et al. 1994), and 
more recent ly non-thermal d et ection s have been 
claimed bv iRephaeh fc Grubeii (1200^ w ith RXTE 
and by iFusco-Femiano et al.l (|1999l I2004D with Bep- 
poSAX, though the latter detection is controversial 
(|Rossetti fc Molendil 120041: iFusco-Femiano et al.l 120071 ). 
Due to the large field of view (FOV) of these non- imaging 
instruments and the simple characterization of the ther- 
mal gas, the source of this emission remains uncertain. 
Even more recently, long (~ 1 Msec) observations with 
INTEGRAL have imaged extended diffuse hard X-ray 
emission from Coma, though it was fou nd to be com- 
plete l y consistent with t hermal emission (IRenaud et al.l 
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Most recently. 



120071: iLutovinov et al.l l2008ir 



Wik et al.r(|2009[) performed a joint 



analysis of spectra from the XMM-Newton EPIC-pn and 
Suzaku HXD-PIN instruments - the most sensitive in- 
struments at soft and hard energies to date - of the Coma 
cluster and were unable to detect IC emission. Instead, 
th ey found an upper lirnit 2.5 times below the detections 
of iRephaeli fc Grubeii (|2002D and IFusco-Femiano et al.l 
()2004[) . However, the narrower FOV of the HXD rela- 
tive to the colhmators of the RXTE PCA/HEXTE and 
Beppo-SAXVDS leaves open the possibility that the spa- 
tial distribution of IC photons is highly extended, and 
therefore much of the flux was missed by the HXD. The 
IC would have to be much broader than the size of the ra- 
dio halo. A uniform IC surface brightness of at least 30' 
in radius from the cluster center is sufficient to reconcile 
these results. Therefore, an imaging analysis at hard X- 



rays is required to confirm this picture; unfortunately, no 
focussing hard (>10 keV) X-ray telescope has yet been 
deployed. In the meantime, it is possible to perform a 
crude imaging analysis with coded ma s k inst ruments, as 
previously discussed bv IRenaud et al.l ([20061 ). 

In this work, we report on the spatial and spectral 
hard X-ray emission from the Coma cluster using the 58- 
month accumulation of the Swift Burst Alert Telescope 
(BAT) all-sky survey. After the first 9 mont hs of the sur- 
vey. C oma was seen to be clearly extended (l Aiello et al.l 
[200l . so an accurate measurement of its flux must ac- 
count for its resolved nature; the standard method of 
extracting fluxes from coded mask instruments assumes 
the underlying source to be point-like. Using models for 
the spatial distribution of thermal and potential non- 
thermal emission, we measure the total, extended flux 
in the 8 energy bands that make up the survey. These 
fluxes are then converted into spectra, which we jointly 
flt with an XMM-Newton EPIC-pn spectrum from a spa- 
tially identical region. In this way, despite poor spatial 
resolution (^ 20'), we are sensitive to any large-scale, 
extended emission above the detection threshold for the 
survey. While the sensitivity of the BAT detector is lower 
than instruments such as the Suzaku HXD-PIN, the sur- 
vey's large exposure time - thanks to a FOV that sees 
1 /8"^ of the sky in a single pointing - gives it a compa- 
rable, if not superior, overall sensitivity. In Section [2l we 
describe the Swift BAT survey in general and the XMM- 
Newton EPIC-pn and BAT observations of the Coma 
cluster specifically. The extraction of spatially extended 
fiuxes from models, along with the specific models them- 
selves, is discussed in Section [3| Spectra constructed 
from these spatial fits are presented in Section [4l along 
with the results of joint fits with the XMM-Newton spec- 
trum. In Section [5l we provide upper limits on spatially 
extended, nonthermal emission, and in Section [6] we dis- 
cuss the implications of our non-detection for the rela- 
tivistic phase of the ICM of the Coma cluster. In the ap- 
pendices, we describe the calibration of the survey such 
that joint fits with XMM-Newton are straightforward, 
and we demonstrate that the BAT instrument intrinsi- 
cally detects extended emission on the scales of interest 
here, though with higher uncertainty than for a point 
source. We assume a fiat cosmology with SI a/ = 0.23 and 
Hq = 72 km/s/Mpc and a luminosity distance to Coma 
of 98.4 Mpc. Unless otherwise stated, all uncertainties 
are given at the 90% confidence level. 

2. OBSERVATIONS 

To achieve the necessary spatial coverage and spectral 
sensitivity, we take advantage of mosaics of the Coma 
cluster constructed from observations by the XMM- 
Newton and Swift satellites. The high sensitivity and 
good spectral and spatial resolution of the XMM-Newton 
EPIC-pn data act as a check on the interpretation of the 
Swift BAT data, allowing the thermal and potentially 
nonthermal emission at hard energies to be accurately 
decoupled. 

2.1. XMM-Newton EPIC-pn Mosaic Observations 

The observations and processed XMM-Ne wton data 
used herein are identical to that presented in lWik et al.l 
|200a) . where a more detailed description can be found. 
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Fig. 1. — XMM-Newton EPIC-pn 2-7.5 keV w avelet-smoothed 
X-ray surface brightness image (.Schuecker et al. 2004) with con- 
tours from the raw Swift BAT 14—20 keV survey image (square 
root spacing: 0.0, 2.625 X lO"*^, 1.05 X lO"'', 2.3625 X 10"^, and 
4.2 X 10~^ counts pix"-*^). Negative contours are not shown 
for clarity, and note that the FWHM of the BAT PSF is 19.'5. The 
(red) box shows the 65f5 X 65.'5 region from which the EPIC-pn 
spectrum is extracted for joint fits. The BAT emission is slightly 
more extended to the W-SW, as would be expected from the higher 
temperature gas in that direction. Note, however, that the outer- 
most contour is consistent with noise and should be ignored. 

The XMM-Newton EPIC-pn mosaic of Coma consists 
of 14 individual pointings, t he first set ( 11 po intings) 
of which were discussed in iBriel et all (|200lD . The 
full 14 pointings consi dered here were first presented in 
iSchuecker eFall (|2004D . For joint fitting with Swift BAT 
spectra, we extract events from a 65.'5 x 65f5 box cen- 
tered on the radio halo, which was originally chosen to 
match the Suzaku Hard X-ray Detector field of view; the 
region is shown a s the outermost contour in Figure 3 of 
iWik et all ()2009[ ). and also as the box in Figure [1] in the 
present paper. The only modific ation o f our analysis pro- 
cedure compared to that in Wik et al.l ([20 09) is that the 
XMM-Newton spectra were not weighted by the spatial 
response of the Suzaku HXD, since they are not being fit 
simultaneously with that instrument. No similar weight- 
ing is needed to comparison to the Swift BAT data since 
the BAT survey covers the entire sky and the vignetting 
of individual pointings is corrected for during processing. 

We also use the temperature map derived from the 
XMM-Newton mosaic and described in I Wik et iail ()2009l ) 
to model the spatial distribution of hard X-ray emis- 
sion, which is detailed in Section [5.1.1l Spectra extracted 
within a 16 X 16 contiguous grid of regions of size 4.'3 x 4.'3 
are fit to single temperature APEC models over the energy 
range 0.5 < < 14 keV. While energies of 9.5 < i? < 14 
keV are generally included during spectral fits, these data 
are primarily used to determine the background and are 
only included as a check of the background subtraction. 
Because of the low count rate at these energies, their 
inclusion does not significantly contribute to or the 
resulting best-fit parameters. 

While data from the Suzaku observation of Coma (OB- 
SID 801097010) are not part of the current analysis, we 
do make use of the 2-7.5 keV spectrum from the chip 
of the X-ray Imaging Spectromete r (XI SO) for calibration 
purposes, as described in Section IAT] We use the same 



spectrum that served to cross-calibrate the data in the 
previous XMM- Newton/ Suzaku analysis (region 10 from 
IWiket al.l[2009l) . 

2.2. The Swift BAT 58-Month Survey 

The Swift mission is primarily to detect and localize 
gamma-ray bursts, which is accomplished with the very 
large FOV (~ 1/8*'' of the sky), coded mask aperture 
Burst Alert Telescope (BAT). As such, the nearly ran- 
dom pointing strategy culminates in an almost uniform 
coverage of the entire sky with an ^ 8 Ms of exposure 
time, made from many ~ 5 minute individual obser- 
vations. Images of the sky are reconstructed by cross- 
correlating the shadow pattern of the randomly coded 
mask in front of the detectors with a detector plane image 
via a fast Fourier transform. The detectors are sensitive 
to hard X-ray/soft gamma-ray photons from 14-195 keV 
in 80 native energy channels. As part of the default sur- 
vey processing, the channels are combined into 8 broader 
energy bands: 14-20 keV (El), 20-24 keV (E2), 24-35 
keV (E3), 35-50 keV (E4), 50-75 keV (E5), 75-100 keV 
(E6), 100-150 keV (E7), and 150-195 keV (E8). The final 
survey is built from the individual sky reconstructions, 
which are summed and resampled onto predetermined 
image planes of 6 facets, each in the Zenith Equal Area 
projection. The detailed processing methodology and 
survey properties for the 58-month BAT all sk y survey 
are ne arly identical to those described in Tuelle r et al.l 
(|2010( ) for the 22-month survey. The only major differ- 
ence is that for the 58-month survey, the gain of each de- 
tector pixel was individually calibrated with an onboard 
radioactive source, which had not been done previously. 
This better accounts for the sensitivity of low gain pix- 
els, effectively increasing the overall sensitivity. Also, the 
sky images are more finely resampled so that the pixels 
near the center of the image projection scale 2f8 instead 
of 5' as with previous versions of the survey. The main 
advantage of this change is to improve the centroiding of 
sources. 

Because the systematic uncertainties in the survey- 
averaged spectrum of the Crab Nebula are smaller than 
the uncertainties in the BAT survey response matrix, 
BAT sur vey fluxes are tied to the Crab fluxes in each 
band (see lTueller et al.ll2010l Sec. 4.5). One drawback to 
this approach is that a source flux is only guaranteed to 
be correctly determined if its spectrum within the energy 
band is identical the Crab's (a power-law with a photon 
index of F ~ 2.1), since the energy response within the 
band may not be uniform. Because we will fit the BAT 
spectrum jointly with the XMM-Newton spectrum, it is 
important that the cross-calibration between the Swift 
BAT and XMM-Newton be accurate. Since the fiux cal- 
ibration of the BAT survey is based on the Crab spec- 
trum, we have made sure that the Swift BAT and XMM- 
Newton agree on the flux and spectrum of the Crab. The 
cross-calibration between Swift BAT and XMM-Newton 
is discussed in detail in Appendix [Xj 

2.2.1. The Coma Cluster 

The 6 facets of the BAT all sky survey are oriented 
in Galactic coordinates with one facet centered on each 
of the Galactic poles and the other 4 centered uniformly 
around the Galactic plane. The Zenith Equal Area pro- 
jection conserves surface brightness but not shapes, so 
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Fig. 2. — The radial profile of the Coma cluster (crosses, black) 
compared to a point source of comparable brightness (squares, red; 
scaled slightly to match Coma's central flux). Each point repre- 
sents an individual pixel. The BAT emission from Coma is clearly 
extended and not axially symmetric, as shown by the the larger 
spread in pixel count rates in its profile compared to the point 
source. 

objects far from the center of the projection can be some- 
what distorted. However, the fortuitous location of the 
Coma cluster near the North Galactic pole, and thus 
the center of its facet, means that any such distortions 
are negligible. Nevertheless, for all parts of the analy- 
sis image pixels are referred to in terms of their Galactic 
coordinates so that any image projection effects are com- 
pletely eliminated. 

Hard X-ray emission is clearly detected in the first 4 
energy bands, up to 50 keV. In Figure [U we present the 
wavelet-smoothed 2-7 keV XMM-Newton EPIC-pn im- 
age of the Coma cluster mosaic overl aid with contours of 
the Crab-normahzed BAT flux (see iTueller et al.ll2010l 
for a description), which shows the hard X-ray emission 
to be elongated in the same East - West direction as the 
softer emission. As first noted bv lAiello et al.l ()2009[ ) in 
the BAT survey, the Coma cluster is partially resolved 
by the BAT, which is explicitly shown in Figure [2] This 
figure compares the radial surface brightness profile of 
Coma with that of a nearby point source with a compa- 
rable flux. Coma is clearly extended. The points plotted 
are individual pixels. The greater width of the distribu- 
tion for Coma indicates that its surface brightness is not 
circularly symmetric. As shown in Figure [U both the 
BAT and XMM-Newton X-ray emission is elongated in a 
ENE-WSW direction. 

While coded masks instruments have some difficulty 
detecting extended emission, this is only true for emission 
extended on very large scales, when the shadow pattern 
of the mask on the detector plane becomes sufficiently 
blurred. Conceptually, as long as at least 1 detector pixel 
is always shadowed by the mask over the entire extent of 
emission, no flux can be mistaken for the unmodulated 
background. The random placement of mask elements 
results in occasional, continuous clumps of over 5 ele- 
ments across, which means that over a couple degrees 
all of the source flux can be recovered since one mask 
element subtends, on average, about 20'. For the BAT, 
we performed simple simulations that suggest that the 
actual scale over which accurate extended fluxes can be 
recovered, in principle, approaches the size of the FOV 
if the spatial distribution of the source is known a priori. 
The practical problem is that flux from each part of an 



extended source adds systematic noise to every other part 
of the source, eventually drowning the signal in noise. 
However, for small extensions this additional uncertainty 
is not overwhelming. In Appendix |B1 we show that the 
Swift BAT provides accurate fluxes for extended sources 
on the scale of interest for the Coma cluster; in particu- 
lar. Fig. [TUl demonstrates that all of the flux of a source 
extending over several degrees is accurately recovered. 

Because nearby pixels in the BAT survey images are 
correlated, it is more difficult to determine the fiux un- 
certainties for extended sources in the BAT than is usu- 
ally the case for X-ray images. For point sources, the 
flux uncertainties can be determined from the fluctu- 
ations in the local background in the BAT survey im- 
ages. For Coma, we calculate the RMS fluctuations in 
the background (<Tbgd) around Coma in an annulus of ra- 
dius 15 < r < 100 pixels (42' < r < 4?67), as is typically 
done for sources in the BAT survey. These values, and 
the flux uncertainties for extended sources like Coma, are 
derived in Appendix [C] 

3. CHARACTERIZATION OF EXTENDED 
EMISSION IN BAT IMAGES 

To extract fluxes for extended sources, we choose to 
test a priori model distributions, as opposed to using a 
method like the "CLEAN" algorithm (Hogbom 1974), 
which reconstructs fiuxes from an unknown underlining 
distribution assuming the PSF shape only. (See Ap- 
pendix IB. 21 for details about this choice.) We repre- 
sent a diffuse source as a collection of point sou rces, 
each of which is convolved by the PSF (Eqn. IBip and 
summed together. Throughout this work, image data 
are fit to these spatial models using the MPFIT algo- 
rithm ([Markwardtl 120091 ) . which performs a Levenberg- 
Marquardt least-squares minimization to converge on 
best-fit parameter values. 

3.1. Model Spatial Distributions of Hard X-ray 
Emission 

Our goal is to detect IC emission from the same elec- 
trons producing the radio halo; however, the electron 
spatial distribution need not follow the radio halo if B 
varies spatially within the ICM (e.g.. iDolag et al.ll2008D . 
Indeed, there is evide nce that the magnetic field in Coma 
declines with radius (jBonafede et al.ir2010t ). which would 
allow for a more extended relativistic electron popula- 
tion visible through IC interactions with CMB photons. 
To accurately search for this potential signal, we need to 
both choose model distributions for this emission and to 
fully model the thermal emission also present in the BAT 
energy bands. 

3.1.1. Tfiermal Models 

Following the success of the XMM-Newton-deiived 
temperature map for explaining the thermal origin of the 
Suzaku RXD-PIN spectrum (Wik et al. 2009), we use the 
same map to predict the spatial distribution of thermal 
emission at hard energies. For each region of the map, 
the flux of the best-flt APEC model is calculated in each 
BAT energy band and treated as a point source at that 
location. Note that because the APEC and MeKaL look- 
up tables are not defined above 50 keV in XSpec, we 
use MeKa to derive the temperature map fluxes in the 
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4 highest energy bands. Then, for each band the 232 
temperature map region "points" are taken together to 
serve as the diffuse model. To compare this or any of 
our diffuse models to the BAT image data, each point 
is assigned the PSF shape with its peak value equal to 
the point flux, and the overlapping PSFs are summed to- 
gether and sampled at the location of the image pixels. 

Thermal emission is detected in the first 4 BAT energy 
bands E1-E4. Since this emission is an extension of the 
X-ray emission which dominates the XMM-Newton im- 
age and its distribution is known, we fix the location of 
the thermal model to the best-fit position of the model 
for the El band data, where the signal-to-noise ratio is 
the highest. 

Note that for the thermal model "fits" to the Swift 
BAT spatial distribution in various bands, only the nor- 
malization of the model in each band is allowed to vary. 
The spatial distribution within each of the bands is com- 
pletely determined from the XMM-Newton data. 

The XMM-Newton data will also contain any nonther- 
mal emission within the XMM-Newton band. Is it rea- 
sonable to use the XMM-Newton temperature map and 
image to determine the spatial distribution of the hard 
X-ray thermal emission? The spectra in the tempera- 
ture map were fit over the energy range 0.5 keV< E < lA 
keV. While the upper limit of this band is fairly hard, the 
low energy limit guarantees that the spectra are domi- 
nated by softer photons. For any sensible nonthermal 
spectrum, the XMM-Newton spectra are dominated by 
thermal emission. In fact, if there is cool, dense gas along 
a given line of sight, the XMM-Newton-hased model may 
actually underestimate the thermal hard X-ray emission. 
In any case, if there is strong nonthermal emission in 
the XMM-Newton spectrum, it will dominate the harder 
BAT energy bands, and will be uncovered in the joint 
spectral fits to the XMM-Newton and BAT spectra (Sec- 
tion m below) . 

3.1.2. Nonthermal Spatial Models 

The Suzaku HXD-PIN upper limit, which is 2.5 times 
below the RXTE and Beppo-SAX detections, only ex- 
cludes those measurements if the IC emission originates 
from a relatively compact region {R < 20'). More ex- 
tended emission of roughly uniform surface brightness, 
however, would be consistent with both the detections 
and upper limit. In most physical models for the IC, it 
is likely that the surface brightness would decline with 
radius; one exception is the KW model discussed below. 
However, there is no single well-established model for this 
decline. Since our object is to test the possibility that 
the difference between our Suzaku HXD-PIN upper limit 
and the RXTE and Beppo-SAX detections is due to the 
extent of the IC emission, we consider the extreme case 
of a uniform surface brightness disk. Thus, we assume 
any nonthermal emission to take the form of a uniform 
brightness, circular disk with a radius R = 25' (R25), 30' 
(R30), 35' ( R35), 40' (R40), 45' ( R 45), o r 60' (R60). 

Recently, iKushnir fc WaxnianI ()2010l ) proposed an- 
other model for the IC emission of the Coma cluster 
which is consistent with both the Suzaku HXD-PIN up- 
per limit and the RXTE and Beppo-SAX detections. In 
this model, the IC hard X-ray emission comes from a 
separate population of electrons from those in the radio 
halo. This new population of electrons are accelerated at 



the virial accretion shock of the cluster at a very large ra- 
dius. These virial shock accelerated electrons lose energy 
quickly, and form a shell of hard X-ray emission, which 
projects on the sky as a ring with nearly uniform surface 
brightness emission in its interior. We will refer to this as 
the KW model. While most of the flux resides in a ring at 
the cluster virial radius, the amount of flux detected by 
an instrument pointed at the cluster center will depend 
sensitively on its FOV. We take a h the model param- 
eter values for Coma directly from IKushnir fc WaxmaiJ 
(2010) when comparing their expected flux to the con- 
straints imposed by the BAT data, though the only rel- 
evant parameter for the spatial distribution is the accre- 
tion shock radius ^200 = 82' 1. The radial distribution 
of flux is simply geometrical in form, assuming an in- 
fi nitely thin shell at th i s radi us; the expression is given 
in IKushnir fc WaxmanI (poTo l. Eqn. 9). 

In reality, it is unlikely the spatial distribution of emis- 
sion would be as regular and axisymmetric as portrayed 
by these models. However, for the spatial extent we con- 
sider relative to the resolution of the BAT, deviations 
from the idealized models will not particularly impact 
our results. 

We do not assume that the center of the nonthermal 
emission distribution from Coma is the same as the cen- 
ter of the thermal emission. Instead, for each of the non- 
thermal models, the model center is placed at 81 different 
positions on a 9 x 9 grid with 2. '5 spacings around the 
centroid of the large-scale thermal emission. The center 
of the nonthermal emission is taken as one of the param- 
eters to be varied in the fits of the spatial and spectral 
distributions below. 

3.2. Spatial Fits to the BAT data 

The spatial models for the thermal emission alone, or 
for the thermal emission plus nonthermal emission, were 
fit to the pixel values in the BAT images in each of the 8 
BAT bands. As noted above (Section r3.1.ip . the center of 
the thermal model was determined by fitting the center 
in the El band. For the thermal model, only the overall 
normalization was allowed to vary. For the nonthermal 
models, the center was varied for a grid of positions (but 
fixed in each individual fit, see Section [5.1.2p . The model 
normalization (i.e., flux) in each of the 8 BAT bands was 
fit independently for the thermal and nonthermal models. 
That is, the spatial fits made no assumptions about the 
spectrum of either type of emission. 

In Figure [21 we present the Swift BAT images of Coma 
in all 8 energy bands (first and third columns) along with 
the thermal model-subtracted residuals for bands E1-E4 
(center column). The spatial distribution of the thermal 
models is represented with the contours in the first col- 
umn. Note that only positively-valued pixels appear in 
the grayscale, which has a square-root scaling, and that 
for each band pure black corresponds to a slightly dif- 
ferent value. Each panel covers an identical 2? 7 x 1?5 
region of the sky. Residuals from the fits are consistent 
with background fluctuations, as can also be seen in the 
radial proflles shown in Figured 

As noted above (Section I3.1.ip . these thermal models 
are solely based on the XMM-Newton temperature map. 
That is, the spatial distribution is completely determined 
from the XMM-Newton data, fit over the energy range 
0.5 keV< E' < 14 keV, and only the total normaliza- 
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Fig. 3. — Images from the 8 energy bands of the Swift BAT survey (first and third columns). The greyscale follows a square root scaling 
from counts/s (white) to > 10~^ counts s~^ pix~^ (black). The contours in the images in the first column show the best-fit thermal 
model for each band, following a square root spacing. For El, the contours range from to 4.2 X 10~^ counts s~^ pix~^ and for E2 and 
E3, the contours range from to 1.7 X 10~* counts pix~^, spanning 5 contours in each case. For E4, the 3 contours represent the 
model surface brightness at 0, 2.5 X 10~^, and 10~® counts s"^ pix~^. The El contours occur at the same levels as shown in Fig.[l] The 
middle column shows the thermal model-subtracted residual images for E1-E4, with the same for greyscale as the data on the left. The 
residuals show that the thermal spatial models are generally well-mapped to the actual data. Note that the local background is also fit for 
and subtracted from the data in the residual images, so the outer fluctuations are not identical to those in the left column. Emission is 
clearly not detected in the 4 highest energy bands (E5-E8), which are shown in the third column. 
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Fig. 4. — Radial profiles of the data and thermal model fits shown in the first column images of Fig. [3] BAT pixels are averaged in annuli 
of 1' width (crosses), as are the model values for each pixel position (histogram); the residuals are plotted below each fit, on the same scale 
as the fit. The structure in the El and E3 residuals could be due to a slightly larger PSF FWHM and/or a true spatial distribution of 
emission that differs slightly from our models; in either case, the effect on the extracted flux would be less than its la error. 
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tion, independently measured in each band, is allowed 
to vary. Ideally, we would use a harder bandpass to de- 
rive the temperature map, especially since it is known 
that such broad band single temperature fits, when per- 
formed on a truly multi-temperature spectrum, produces 
biased results (Nevalainen ct al. 2010). In practice, how- 
ever, the lower count rate at energies above 2 keV means 
that temperatures cannot be determined in regions with 
lower temperature gas or lower flux. We also extract such 
a temperature map, using the 2-12 keV band, and find 
systematically higher temperatures, < 10 % larger than 
the temperatures in the map we use here. Taking this 
map to represent the input spatial distribution, we ob- 
tain statistically worse fits to the data, especially in the 
El band. Because cooler, less bright gas in Coma's out- 
skirts is not included in this temperature map, though it 
is present in the BAT data, the emission at larger radii 
cannot be properly fit. In any case, the fiuxes we derive 
in each BAT band do not significantly depend on the ab- 
solute temperatures in the map, since only relative tem- 
perature differences across the cluster matter, as the total 
normalization is independently fit for in each band. Be- 
cause all the temperatures in the map are increased when 
a higher energy bandpass is used, the fiux in each BAT 
band remains essentially the same regardless of which 
temperature map is adopted. We therefore use the map 
determined from the 0.5 keV< E < lA keV band in or- 
der to include the cooler gas at larger radius. Regardless, 
this thermal model provides an adequate description of 
both the spatial and spectral (discussed in Section 2] be- 
low) properties of the detected emission in the BAT data. 
This success justifies our approach and confirms that ex- 
tended emission is detected with the same efficiency as 
that from the cluster center. 

The good fit of the XMM-Newton-ha.sed thermal model 
for the emission in BAT bands E1-E4 and the lack of 
obvious excess emission in the harder E5-E8 bands sug- 
gests that nonthermal emission is not very strong or ex- 
tended. The nonthermal model with the mos t extended 
emission is the KW model (Kushnir fc Waxm an 2010), 
in which the IC hard X-ray emission comes from a thin 
shell at a very large radius. Following the methodology 
in Section [B.21 we simulate the expected combined ther- 
mal and nonthermal fiux distribution for this model, and 
compare to the actual data. The results for the El BAT 
band are shown in the upper panel of Figure [S] Clearly, 
we do not detect the nonthermal emission expected by 
the KW model. The lower panel shows the thermal plus 
KW model compared to a simulation of the BAT data as- 
suming the distribution actually followed this model. It 
is clear that the statistics in the BAT data would allow us 
to detect the nonthermal emission from the KW model, 
were it present. A more quantitative limit is derived in 
Section [5l 

4. SPECTRAL FITS 

To search for a nonthermal component in the X-ray 
spectrum of the Coma cluster, spectral models are fit to 
the data. The fits were done for the Swift BAT spec- 
trum alone, or simultaneously with the XMM-Newton 
spectrum of the cluster. The BAT spectra were binned 
into 8 spectral channels, given by the 8 standard BAT 
bands E1-E8. To determine the spectrum in each chan- 
nel, the total measured raw photon fluxes in each band 




— 



20 40 60 80 100 

Radius (arcmin) 

Fig. 5. — The top panel shows the BAT El profile and spatial fit 
for a model with both thermal emission and nonthermal emission 
following the KW model (histogram, red). The best-fitted value 
of the normalization of the KW component is actually negative. 
The dash-dot line represents the background level and the dashed 
line (blue) s hows the predicted spatial distribution of flux for the 
KW model IIKushnir L Waxm an''20ro) from 14-20 keV. The inset 
expands the scale of the y-axis above it to highlight the difference 
between the data and the expected flux. In the bottom panel, we 
perform the same fit to simulated BAT data based on the thermal 
plus KW model, including shot and systematic noise comparable 
to that present in the actual data. This shows that the BAT would 
have easily detected a nonthermal component with the spatial dis- 
tribution given by the KW model and the predicted fiux. 

were converted into "true" photon fluxes using the cal- 
ibration determined from the BAT Crab spectrum; we 
also create a generic redistribution matrix to better rep- 
resent models with spectral shapes that differ from the 
Crab (see Section lA.il) . The uncertainties in each chan- 
nel were determined from the flux uncertainty for a point 
source Cbgd (see Section [C.l|) . and then corrected for the 
effect s of source extent as described in Section[C21 Equa- 
tion (|C3p . and Table [3l T he final uncertainty is given by 
f diffuse in Equation (jC3[) . Additionally, a problem with 
the implementation of the APEC, which we use as our ther- 
mal description for spectral fitting, and MeKaL models in 
XSpec is that the look-up tables are undefined above 50 
keV. Therefore, for the 4 energy channels above 50 keV, 
we substitute MeKa for APEC. This should have no signifi- 
cant effect given the small thermal fiux at these energies 
relative to the errors. 

For the BAT-only spectra, in which only the nonther- 
mal component of the spatial fit is used to build the 
spectra, a single power law model is sufficient to mea- 
sure the nonthermal flux. We also do joint fits of the 
Swift BAT (containing both the thermal and nonther- 
mal spatial components) and XMM-Newton spectra. In 
these fits, the excellent statistics at low energies in the 
XMM-Newton spectra very strongly constrains the ther- 
mal emission. However, the XMM-Newton mosaic covers 
a smaller area compared to either the Swift BAT or most 
of the nonthermal spatial models. (The XMM-Newton 
extraction region used for this spectral analysis is indi- 
cated by the square in Figure [TJ) Thus, in these fits the 
models applied to the XMM-Newton spectra are reduced 
by the fraction of the emission in our XMM-Newton spec- 
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tral extraction region. 

4.1. Joint XMM-Newton EPIC-pn - Swift BAT 
Thermal Emission Fit 

We first consider a purely thermal model for the X-ray 
emission in Coma, and fit the Swift BAT and XMM- 
Newton spectra simultaneously. The spatial distribution 
of the emission was assumed to follow the thermal model 
(Section I3.1.ip as determined from the XMM-Newton 
data. For these and all following fits, we restrict our- 
selves to energies above 2 keV. This choice ensures that 
the thermal component is not strongly biased to lower 
temperatures by cooler gas in these global, multi-phase 
spectra (e.g.. iNevalainen et al.ll2010l) . The resulting sin- 
gle temperature fit is presented in Figure [6l and the 
parameters are given in the first row of Table [TJ The 
quality of the fit is quite good, indicating that a single 
component description for the temperature structure is 
sufficient and that the spectra have been reasonably well 
cross-calibrated. A slight ascending trend in the E1-E3 
residuals exists, however, which is primarily due to a 
slightly lower than expected El flux. While not particu- 
larly significant, it is worth mentioning several potential 
causes for the trend. The most straightforward explana- 
tion is that the calibration is slightly wrong. We presume 
the true Crab spectrum to be a simple power law across 
the entire 2-200 keV interval, but if the spectrum actu- 
ally steepens aroun d E ^ 10 keV as is likely the case 
([Kirsch et al."2005V the higher energy bands will be as- 
signed progressively higher fiux conversion factors (col- 
umn 3 in Table [3]) ; basically, the flux in an energy band 
will be more and more overestimated for bands at higher 
and higher energies. Also, because emission is more sig- 
nificantly detected in the lower energy bands, a small 
change in the overall value of the XMM-Newton and 
Swift cross- normalization factor - such that the Swift flux 
would be raised - could reduce the spread in residuals. 
From a more physical perspective, a single temperature 
model is not entirely appropriate; in multi-temperature 
model flts, the trend in residuals is not as strong. In any 
case, an adjustment to the cross-calibration of 3-5% is 
sufficient to account for the trend, which is well within 
our assumed 90% confidence interval uncertainty of 10%. 

The best-fit temperature is kT = 8.24 ± 0.12 (stat) ± 
0.15 (sys) keV. The systematic term in the error is based 
on varying the cross-normalization factor by 10%; the 
origin of this percentage is discussed in Section [S] be- 
low. This global temperature matches e xtremely well 
with p revious measurements. For example. iHughes et al.l 
(|1993D found kT — 8.21 keV with a spectrum span- 
ning ~ 2 < £' < 11 keV from the Ginga satellite, 
which had a simflar F OV (collimato r with 1-2° FWHM) 
to our aperture, and Arnaud et al.l ()2001f ) found kT = 
8.25 keV over a smaller FOV and lower energy range 
(0.3 < < 10 keV) with the XMM-Newton EPIC-MOS 
instruments. Including lower energy photons tends to 
lower single-temperatu re fits to multi-component spectra 
(jCavagnolo et al.ll2008l ) , and having a smaller FOV tends 
to emphasize the hotter central temperature of kT ~ 9 
keV in Coma, which both explains why the Ginga and 
XMM-Newton temperatures agree and why t he XMM- 
Newton- Suzaku analysis of iWik et aTl (|2009[ ) found a 
slightly higher kT = 8.45 keV. While we find good agree- 
ment with other measurements, note that our tempera- 




Energy (keV) 



Fig. 6. — Single temperature (APEC, red histogram) fit to the 
XMM-Newton EPIC-pn and Swift BAT spectra. The BAT spec- 
tra shown are reduced to the fraction which occurs in the XMM- 
Newton spectral extraction region. The BAT spectra were con- 
structed assuming the spatial distribution predicted by the XMM- 
Newton temperature map. The origin of the low El (14-20 
keV) flux is discussed in the text. A single temperature model 
(kT = 8.24 keV) is sufficient to describe the 2-200 keV emission 
from the central square degree region of the Coma cluster. 

ture, along with other temperatures derived with XMM- 
Newton, could be systematically cooler by a few tenths 
of a keV than temperatures obtained with other obser- 
vatories, given the steeper XMM-Newton Crab spectrum 
and the implications for its instrumental response (see 
Section IATT|) . 

We also tried a two-temperature thermal model for 
the XMM-Newton and BAT data. However, the two 
temperatures and normalizations could not be individ- 
ually constrained by the data (Table [T]); the tempera- 
tures/normalizations listed in the table result when the 
two-temperature model is fit for with initial temperatures 
of kTi = 6 keV and kT2 = 10 keV. The two-temperature 
fit is not a significant improvement on a single tempera- 
ture flt. 

While the average spectrum in the square degree re- 
gion around Coma is adequately described with one or 
two temperatures for the gas, in fact the temperature 
distribution is quite non-isothermal. We account for 
temperature variations in the spatial models used to ex- 
tract fluxes from the BAT images by extrapolat i ng the 
XMM-Newton temperature map from iWik et al.l ()2009f ) 
to higher energies. This map can also be converted into 
a spectral model (labeled "Tmap" in Table [1]) and com- 
pared to the joint spectrum. By including the spatial 
information of the XMM-Newton data in the spectral 
model, we can better account for the thermal contribu- 
tion in the BAT energy bands. The quality of this fit 
(allowing the normalization, but not the shape, of the 
model to vary) is reported in Table [1] In principle, the 
"Tinap" spectral model should perfectly represent the to- 
tal XMM-Newton spectrum, but due to incomplete cover- 
age of the temperature map with the total XMM-Newton 
region and the imperfect determination of the individual 
temperatures, the value - which is primarily driven by 
the higher quality XMM-Newton data - is larger than for 
the other fits, in which the model shape is free to vary 
and can account for these minor differences. Even so, 
the BAT data are slightly better described by this model 
than by any of the other spectral models presented. Al- 
though this result is perhaps expected, given that the 



10 



Wik et al. 



TABLE 1 

Joint Fits to XMM-Newton and Swift Spectra 



Spatial 




SpcctrSyl 


kT 


Norm. 


r or kT 


Norm . ^ 


/dof 


Model 




Model 


(keV) 


(cm^^) 


b 






Thermal Region 


Single T 


8.28 ± 0.13 


0.373 ± 0.002 






1576.79/1544 


Thermal Region 




7.8 


0.25 


9.4 


0.12 


1575.77/1542 


Thermal Region 


Tniap 










1590.01/1545 


Thermal Region 


T+IC<= 


8.27 ±0.13 


0.373 ± 0.003 


7.2 


<0.51 


1576.79/1543 


Thermal Region 


T+IC* 


8.27 ±0.13 


0.373 ± 0.003 


2.0 


<0.00148 


1577.21/1543 


Thermal & 


L KW 


Single T 


8.30 ±0.13 


0.373 ± 0.002 






1570.32/1544 


Thermal & 


L KW 




7.8 


0.25 


9.7 


0.12 


1568.63/1542 


Thermal & 


c KW 


Tmap 










1584.59/1545 


Thermal & 


c KW 


T+IC^ 


8.30 ±0.14 


0.372 ± 0.002 


-1.5 


<1.9 X 10-s 


1570.23/1543 


Thermal & 


L KW 


T+IC^ 


8.30 ±0.14 


0.373 ± 0.003 


2.0 


<0.00082 


1570.62/1543 



^ Normalization of the APEC thermal spectrum, which is given by {10~"'"^/[47r(l+2:)^D^]} J nemi dV, 
where z is the redshift, is the angular diameter distance, rie is the electron density, is the 
ionized hydrogen density, and V is the volume of the cluster. 
Value is T for the T+IC model and kT (in keV) for the 2T model. 

Value is the normalization of the power-law component for the T+IC model, which is the photon 
flux at a photon energy of 1 keV in units of photons cm~-^ keV~^. For the 2T model, the value 
is the normalization of the second APEC thermal model in units of cm~^. 

Parameters unconstrained. 

Value of r is fixed when deriving errors. 
' Value of r fixed based on radio spectrum. 

"Tmap" spatial distribution is used to measure the BAT 
fluxes, it does indicate that the method is self-consistent. 

In all of the thermal models investigated, no evidence 
for a strong high-temperature component is hinted at by 
the BAT data. 

4.2. Nonthermal Spectra 
4.2.1. Nonthermal Emission from the Cluster Center 

To search for evidence of more centrally located non- 
thermal emission, an IC component was first fit to purely 
thermally-derived spectra - i.e., spectra created from fits 
to the BAT data using only the thermal spatial model 
- which arc reported in rows 4 and 5 of Table [TJ As in 
Section |4T1 all fits are to the joint XMM-Newton- Swift 
spectrum extracted from the square region in Figure [TJ 
Not surprisingly, the thermal model parameters are al- 
most identical to the fits without the IC component, and 
no significant IC emission is present. The good single 
temperature fit to the Swift BAT and XMM-Newton data 
already suggests that the nonthermal contribution is not 
very significant. If the photon spectral index F of the 
nonthermal component is allowed to vary, it is uncon- 
strained and the best-fitted value is unphysically steep, 
and in any case, only an upper limit can be placed on the 
nonthermal flux (Table [1} . If we assume a photon index 
of 2.0 for the nonthermal component based on the radio 
data {Section[5]below), the 90% upper limit on the 20-80 
keV flux is 1.24 x lO"""^^ ergs/s/cm^, which corresponds 
to < 0.8% of the total flux in the range 2 < < 200 
keV. 

4.2.2. Extended Nonthermal Emission 

Our goal is to search for extended nonthermal emis- 
sion, which is measured with spatial model fits to the 
unlimited FOV BAT survey images (Section 13. 2p . Be- 
cause the thermal and nonthermal spatial model nor- 
malizations are individually and simultaneously allowed 
to vary to best match the total flux in the BAT images, 
spectra can be created from the sum of both components 



or separately, and also within any aperture. For each of 
the 81 grid positions relative to the center of the cluster 
at which each nonthermal model was fitted for, two spec- 
tra are produced. One consists of the total, unvignetted 
fiux of the nonthermal component only, and the other 
includes the sum of both the thermal and nonthermal 
emission inside the XMM-Newton extraction region. The 
second type of spectrum has the advantage that it can 
be jointly fit with the EPIC-pn spectrum, which in the 
case of a non-detection provides a tighter constraint on 
the flux of nonthermal emission than the flrst type, since 
the nonthermal component must be consistent with the 
higher quality, lower energy XMM-Newton data as well. 
In none of these cases, for either type of spectrum, is a 
nonthermal component detected with > 2(t confldence. 
Wc therefore conclude that, while the Swift BAT instru- 
ment is certainly sensitive to extended emission, none of 
a nonthermal origin is observed in the current version of 
the survey. As an example, the fit parameters for var- 
ious spectral models are shown in Table [T] for the KW 
spatial model nominally positioned (i.e. centered on the 
large-scale XMM-Newton emission). 

In Figure [71 the nonthermal model with the most sig- 
nificant IC component is shown, assuming a fixed photon 
index F = 2 for the spectral fits of each nonthermal spa- 
tial distribution tried. Note that the model in this figure 
represents the upper limit for a nonthermal component, 
not its best-fit value, and that the cross-normalization 
has been adjusted by 10% in the direction that favors a 
nonthermal signal. The BAT spectra in Figures [5| and [7| 
are quite similar, indicating that even in the most favor- 
able case the data reject a significant IC contribution to 
the spectrum of the Coma cluster, extended or otherwise. 

5. UPPER LIMITS TO DIFFUSE, NONTHERMAL 
EMISSION 

To ensure appropriate upper limits are derived, we 
must determine and include any important systematic 
uncertainties in our results. Typically, instrumental 
and/or cosmic backgrounds can be a serious concern and 
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Fig. 7. — Thermal (red histogram) plus nonthermal (blue his- 
togram) model fit to the Swift and XMM-Newton data for the 
r = 2 power law nonthermal model corresponding to the 90% up- 
per limit, including the systematic uncertainties as described in the 
text. This example, which is the model with most significant non- 
thermal flux, is for a 25' radius, uniform surface brightness disk of 
nonthermal emission with a position offset from the center of the 
large-scale thermal emission by —2'.5 and 10' in I and b, respec- 
tively. As is true of all the upper limits from the joint spectra, the 
nonthermal component does not exceed ~ 1% of the low energy 
part of the XMM-Newton spectrum, nor does it compete with the 
thermal emission until energies E > 50 keV. 



must be carefully treated. At lower energies where the 
thermal emission is bright, the background is not compa- 
rable to cluster emission until an energy of '-^ 7 keV, so 
even a background uncertainty of a few percent does not 
significantly impact the XMM-Newton EPIC-pn spec- 
trum. This background includes both the non- X-ray and 
cosmic backgrounds. Given the size of the mosaic, the 
uncertainty in the overall EPIC-pn background used here 
(-2-12 ke'v) is 2.4% (Read & Ponnian 2003). Increasing 
the XMM-Newton background favors larger nonthermal 
fluxes since this reduces the temperature slightly, so the 
XMM-Newton background is raised by 2.4% during up- 
per limit derivations. Point sources in the XMM-Newton 
mosaic account for only — 1% of the emission and have 
a spectrum that as a w hole does n ot vary significantly 
from a F ~ 2 power law (IWik et al.l[200a) . so we do not 
model their contribution to the XMM-Newton spectrum. 
While their flux may artificially enhance a nonthermal 
signal, ignoring them will only result in slightly more 
conservative upper limits. For the BAT survey data, the 
background is automatically removed as part of the im- 
age reconstruction procedure, and systematic variations 
are encoded as fluctuations in empty sky regions, which 
is already included in the error budget. 

The more significant systematic uncertainty is in the 
determination of the cross-calibration between Swift and 
XMM-Newton. Ideally, there should be no uncertainty 
since we based the BAT calibration on the XMM-Newton 
data and the Suzaku XISO Crab spectrum. However, the 
slope of the calibration (i.e. the assumed photon index 
of the Crab) is less certain. The total error, statistical 
and systematic, of the p hoton index in XMM-Newton 
EPIC Crab fits is ±0.05 (jKirsch et all [2005). so we ad- 
just our calibration to make the canonical Crab spectrum 
flatter by 0.05, which acts to increase the BAT fluxes 
(~ 10% for El, - 20% for E8) and flatten the BAT spec- 
tra, thus enhancing nonthermal fluxes. This approach to 
the systematic uncertainty is also conservative, as it is 



know n that the C rab spectrum steepens above 10 keV 
(e.g., iKirsch et al.] [2005). 

Since we have no clear detection of nonthermal emis- 
sion, we must decide on its photon index from other ar- 
guments. The natural choice is to use the spect ral in- 
dex o f the radio halo, or F = 1.5-2.5 ( [Giovannini et al.l 
I1993D . though the lower range of the X-ray regime ex- 
plored here corresponds to lower energy electrons where 
the emission may have a flatter spectrum. Also, both pre- 
vious detections using RXTE and Beppo-SAX data found 
F ~ 2, and the model of Kushnir & Waxman (2010) pre- 
dicts this photon index. Therefore, we fix the nonther- 
mal power law index to F = 2, primarily because we are 
most interested in directly comparing our upper limits 
with these previous detections and model predictions. If 
the spectrum of nonthermal emission is in fact flatter, 
our upper limits will be low by some amount since the 
BAT errors are large and the XMM-Newton data will 
have less leverage on the fits. However, the high energy 
flux will not increase d ramatically: as illustrated in Table 
3 of IWik et al.l (|2009t) . the 20-80 keV flux rises by a fac- 
tor of 2 from F = 2tor = 1.5, and trials show the same 
behavior for the nonthermal component in this work. 

For the above systematic uncertainties, we flnd the 
90% confldence upper limits to nonthermal emission 
for both the nonthermal-only BAT and for the joint 
EPIC-pn/BAT spectra. The thermal component in 
the latter case is simultaneously fit with the normal- 
ization of the nonthermal component. We present 
each individual limit in Figure [5] along with the 
vignetting corrected flux es / upper limit from RXTE 
(iRephaeli fc Gruberl 120021 upper cross, green), Beppo- 
SAX (iFusco-Femiano et al.ll2004 lower cross, red), and 
Suzaku C Wik et al.l l2009'. upper limit, blue). The colli- 
mator responses for these instruments are convolved with 
the model flux distributions to give these values or limits. 
The Suzaku HXD-PIN instrument is a square collimator 
wit h spatial sensitiv ity of the form given in Equation 3 
of (jWik et al.l 120091) ■ and the RXTE PCA/HEXTE and 
Beppo-SAX PDS instruments are hexagonal collimators 
with triangular approximation FWHM of 1° and 1?3, re- 
spectively. We approximate the nearly axisymmetric re- 
sponse with a 4th-order polynomial of the form: 



i?hex(0) = 1.00 - 1.36 (^) +0.46 (jf;^)' 



-0.58 



0.48 



(3) 



where i?hcx(^) is the fraction of emission visible to the 
instrument at off-axis angle and 0max is where emission 
is no longer detected. For each spatial model, our upper 
limits are ordered in Galactic coordinates from the low- 
est values of I and b in our grid, incrementing I for all 
positions with that latitude before incrementing 6, with I 
reset to the minimum value. It is this ordering that pro- 
duces the pattern evident in the limits. The limits for 
the joint spectral fits are given in the top panel, while 
the nonthermal-only spectral limits are provided in the 
bottom panel. 

Surprisingly, the upper limits derived from extended 
spatial models and the BAT data alone (bottom panel) 
are comparable in sensitivity to the previous detec- 
tions/limit. Larger models are generally less constrained, 
due to the greater uncertainty in their estimated flux 
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Fig. 8. — Upper li mits (sm all arrows) for each nonthermal spatial model relative to the RXTE f Rephaeli fc Gruberll2003 . upper cross, 
green), Beppo-SAX IjFusco-F emiano et al. 2004, lower cross, red), and Suzaku (Wik ct al. 2009, long arrow, blue) detections/upper limit. 
The previous flux detections [(f .5it0.5) X 10^^^ ergs/cm^/s] and upper limit (6x 10~^^ ergs/cm^/s) are corrected to account for the fraction 
of emission missed due to vignetting by the collimator response functions (see text). In the top panel, upper limits are calculated from 
the simultaneous joint fits to the XMM-Newton and Swift spectra, and all nonthermal fluxes reported are from inside the XMM-Newton 
extraction region (the square in Fig. [ij. In the bottom panel, upper limits are derived from BAT spectra created from the nonthermal 
component of spatial fits only, and the fluxes represent the total emission of the spatial model. Based on the results presented in the top 
panel, we conclude that extended IC emission cannot reconcile the discrepancy between the Suzaku and RXTE/ Beppo-SAX observations. 



(Equation IC3p . though local fluctuations have a greater 
impact on smaller models, increasing the spread with 
position. Stronger constraints are obtained when lower 
energy emission is simultaneously considered (top panel), 
and a similar range in upper limits is found for each spa- 
tial model distribution. This result is not surprising, as 
each model contributes roughly the same amount of flux 
inside the XMM-Newton extraction region, since most 
of them are extended beyond this region. Interestingly, 
the KW model provides the limits most consistent with 
the Beppo-SAX detection, which follows from the large 
PDS FOV - it would observe a higher proportion of the 
brighter ring emission - and the larger errors resulting 
from the KW model's size. However, in all realizations 



of the joint fit case, our 90% limits lie below the 90% 
interval of the previous detections/limit. Thus, all of 
the previous detections are excluded for any of the spa- 
tial models when one fully accounts for the differences 
in spatial sensitivity between instruments. Our upper 
limits, for the nominal case where the nonthermal dis- 
tribution is centered on the large-scale thermal emission, 
are compared to these previous measurements in Tabled 

6. IMPLICATIONS AND DISCUSSION 

By taking advantage of the crude imaging capabilities 
of the Swift BAT instrument and the impressive sensitiv- 
ity of the 58-month all sky survey, we are able to con- 
strain the amount of nonthermal, hard X-ray emission - 
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TABLE 2 

Flux Upper, Limits (20-80 keV) for Nominal NT Position'' 



Model 


Joint BAT limit 


Beppo-SAX 


RXTE 


Suzaku 




10-12 cgs 


10-12 cgs 


10-12 cgs 


10-12 cgs 


KW 


< 3.86 


16.8 ±5.6 


8.3 ±2.8 


< 6.2 


R60 


< 2.16 


19.3 ±6.4 


13.4 ±4.5 


< 25.1 


R45 


< 2.34 


23.0 ± 7.7 


18.4 ±6.1 


< 25.4 


R40 


< 2.53 


25.3 ±8.4 


21.0 ± 7.0 


< 24.7 


R35 


< 2.95 


26.1 ± 8.7 


22.4 ±7.5 


< 21.7 


R30 


< 3.22 


23.7 ± 7.9 


21.0 ± 7.0 


< 16.8 


R25 


< 3.48 


21.7 ± 7.2 


19.7 ± 6.6 


< 13.3 



^ Units for flux are 10 12 erg cm 2 g l 

extended or otherwise - from the Coma cluster. We find 
no evidence for an extended, hard excess tha t could rec- 
oncil e recent detections from RXTE (Rephaeli fc Grubeii 
I2002D and Beppo-SAX (iFusco-Femiano et al.ll2004D with 
the upper hmit from Suzaku ( Wik et al.l I2OO90 : note, 
however, these det ections would still b e in co nflict with 
the upper hmit of iRossetti fc Molendil ()2004D . Generic, 
uniform surface bri ghtness disks, along with a recently 
proposed IC model (jKushnir fc Waxm"anll2010f) . were fit 
to BAT survey images, converted to spectra, and investi- 
gated for signs of a nonthermal component. For each spa- 
tial model, we compute upper limits on a grid of positions 
and compare them to previous measurements, being care- 
ful to convert detected fluxes into intrinsic source fluxes, 
given a particular spatial distribution, by accounting for 
the collimator vignetting functions. These are direct 
comparisons, in the sense that the instrumental response 
of all detectors involved have been fully co nsidered, and 
as such we, like IRossetti fc Molcndi (2004) , cannot con- 
firm the_claimMd^tection of Rephaeli fc Gruber (2002) 
and iFusco-Femiano et all (jl004). 

Their observed hard excesses could have had other 
reasonable sources, if not diffuse IC emission from the 
nonthermal phase of the ICM. A common difficulty is 
an accurate determination of both the cosmic and non- 
X-ray background, the treatment of which is the pri- 
mary difference between["Fiisco-Femiano et al.l (|2004) and 
IRossetti fc Molendil (j2004D . Another possibihty is the 
variable nature of nearby point sources, most notably 
the AGN W Comae, which was once quite bright but 
has been fading for many years. The concurrence of the 
RXTE and Beppo-SAX observations could have simply 
been unlucky and caught W Comae (or another source) 
in a bright state. A somewhat more subtle, and perhaps 
more likely, explanation concerns the multi-temperature 
nature of Coma's ICM. Small amounts of hot gas could 
dominate the high energy emission, so the extrapola- 
tion of an average temperature determined from lower 
energy data may not be an adequate description of the 
thermal contribution to the high energy flux. The ef- 
fect of the multi-temperature gas in Coma is evident in 
the SW extension at hard energies observed by INTE- 
GRAL (|Renaud et al.l[2006l : lEckert et al.][200l and con- 
firmed here; higher tem peratures seen at this location 
in the temperature map (|Wik et al.ll2009t l are sufficient 
to explain the change in morphology, which points to 
the increased significance of this gas at higher energies. 
Even so, single temperature spec tral fits do not produce 
IC detections in this study or in lWik et ahl ([2009. ) . The 



explanation may simply rest in a slight mischaracteriza- 
tio n of the hard energy e missio n weighted temperature; 
in IFusco-Femiano et all (|2004l ). the FOV of the lower 
energy HPGSPC instrument does not quite match the 
higher energy PDS , and the tem p eratu re of 7.67 ±0.1 
keV found in Rep haeli fc Grubeii (|2002f ) is significantly 
below that allowed by the XMM-Newton data. 

Given radio synchrotron emission and an upper limit 
on the X-ray IC flux, a lower limit on the average ICM 
magnetic field can be estimated, as described by equa- 
tion (13) in Wik et al. (2009) and the accompanying text. 
A diffuse radio flux of 640 mJ y at 1.4 GHz is d etected 
out to a radius of 40' in fPeiss et all (|1997D . For 
comparison, we will use the upper limit of 2.7 x 10^^^ 
erg s-i cm-2 {20 < E < 80 keV) from the R = 40' 
disk model. These val u es im ply B > 0.25 /iC, an in- 
crease from iWik et al.l (|2009D but still well below the 
equiparti tion value of Bp„ = 0. 5 /iC for the Coma ra- 
dio halo (|Giovannini et all 1199 3^. A slightly lower limit 
of i? > 0.2 ^G results if a more conservative IC upper 
limit of 4.2 x 10~^^ erg s~^ cm-^ is used, which consid- 
ers the limits from all spatial models tested. Regardless, 
these limits on B fall well below line of sight estimates 
of several nG from Faraday r otation measure (RM) ob- 
servations (jFeretti et al.lll995l) . though due to geometric 
effects these measurement s may not repre sent the aver- 
age cluster magnetic field (i Petrosianll2001l ) . 

However, the global field may be recovered by com- 
bining many RM measurements along different lines- 
of-sight through th e I CM with nume r ical s imulations 
(|Murgia et al.l l2004[) . IBonafede et"all (j2010D have ap- 
plied this method to the Coma cluster, deriving a radial 
profile where the energy density of the magnetic field falls 
roughly in proportion with the energy density of thermal 
gas and with a central field strength of Bq ~ 4.7 /iC. 
Combining this model of B[r) with an approximate rep- 
resentation of the radial density profile of synchrotron 
emission, implied by a ro ugh /3-model fit to the point 
source subtracted image of iDeiss et al.l (|1997l ) (tc = 18', 
/3 = 1, and /q = 1.23 mJy arcmin-^), directly leads to 
a prediction of the expected IC surface brightness as a 
function of radius. Our illustrative - due to the large un- 
certainties in all parameters assumed in this exercise - IC 
surface brightness distribution is flat out to ~ 30' with a 
20-80 keV flux of ~ 8 x 10^^'' erg s~^ cm-^ arcmin-^, at 
which point it nearly linearly drops toward zero, though 
not reaching it, around a radius of 90'. This surface 
brightness is about an order of magnitude below that 
implied by our upper limits, providing a possible expla- 
nation for why we are unable to detect an IC signature. 
An order of magnitude lower hard X-ray flux for Coma 
is also predicted bv .Brunetti fc La zarian (2010, see their 
Fig. 5) , who have developed the most comprehensive pic- 
ture yet of radio halo generation by MHD turbulence. 
On the other hand, larger IC fluxes would be expected 
if the radio synchrotron emission falls off more gradually 
than modeled here, since a flatter radial profile would 
suggest a higher relativistic electron density given the 
falling magnetic field with cluster radius. More accurate 
maps of Coma's radio halo, preferably at lower frequen- 
cies where the radio electrons correspond more closely to 
the IC-emitting electrons, will clarify this issue. 

Ultimately, a true detection of IC emission from Coma 
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will have to wait for upcoming missions with focussing 
hard X-ray telescopes, namely NuSTAI^ and Astro-H 
Fl . For NuSTAR to achieve a sensitivity comparable to 
our upper limits, a singl e pointed observatio n of at least 
100 ks will be required (Madse n et al.l[2"009f ). However, 
the much finer spatial resolution will remove the uncer- 
tainty associated with bright background AGN and al- 
low multiple spatially-resolved joint fits. Assuming the 
hottest gas, which produces the largest amount of ther- 
mal emission at hard energies, is localized, then these 
regions can be identified and avoided in order to de- 
tect a lower surface brightness, but more uniform, IC 
component. Similarly, if the IC emission is more local- 
ized, it will be easier to identify with spatially-resolved 
joint fits between XMM-Newton and NuSTAR or Astro- 
H spectra, as has bee n done with Chandra data alone 
(jMillion fc Alle n '20091. The unambi guous detection of 
IC emission associated with radio halos and relics is cru- 



cial to determining the energy content in the relativistic 
phase of the ICM and how significant of an influence this 
phase has on the dynamics and structure of the thermal 
gas in clusters. 

W e thank W . Reich who kindly provided us with 
the iDeiss et al.1 ()1997[ ) radio image, C. B. Mark- 
wardt who explained to us (and wrote) many of the 
BAT software analysis routines used in this work, 
and the referee whose comments helped clarify the 
discussion of several key points. DRW was sup- 
ported by a University of Virginia GSAS Dissertation 
Year Fellowship and a Virginia Space Grant Consor- 
tium Fellowship. DRW and CLS were supported in 
part by NASA through Suzaku grants NNX08AZ99G, 
NNX09AH25G, and NNX09AH74G, and XMM-Newton 
grants NNX08AZ34G and NNX08AW83G. Basic re- 
search in radio astronomy at the NRL is supported by 
6.1 Base funding. 



APPENDIX 

A. XMM-Newton EPlC-PN-Swift BAT CROSS-CALIBRATION 

Because the systematic uncertainties in the survey-averaged spectrum of the Crab Nebula are smaller than the 
uncertainties in the BAT survey response matrix, BAT survey fluxes are tied to the Crab fluxes in each band since 
the systematic uncertainties in the surve y-averaged spectrum of the Crab Nebula are smaller than the uncertainties 
in the BAT survey response matrix (see iTueller et all 120101 Sec. 4.5). This method also requires that the intrinsic 
Crab spectrum be defined since its exact spectrum remains somewhat uncertain (see Sec. IA.2l below'). particularly at 
higher X-ray energies. In practice, though, we are less concerned with an accurate absolute calibration for the BAT 
than we are with, in this case, an accurate calibration relative to our XMM-Newton EPIC-pn spectrum. Therefore, 
instead of prescribing a canonical Crab spectrum as close to the true spectrum as it has been measured thus far, we 
need to set it to the Crab spectrum as measured by the XMM-Newton EPIC-pn instrument over the energy range 
we consider. Otherwise, systematic calibration errors between the instruments could significantly affect our result, 
since our goal is to detect excess radiation at hard energies due to nonthermal emission. Errors leading to steeper 
(flatter) XMM-Newton spectra and flatter (steeper) Swift spectra, for example, will reduce (increase) the thermal 
contribution at higher energies and similarly enhance (suppress) a nonthermal signal. In other words, any systematic 
miss-calibrations are mimicked in the BAT calibration so that thermal and nonthermal models can be simply applied 
during joint fits of the data. 

A.l. The Spectrum of the Crab According to XMM-Newton 

Because of its high X-ray flux, simple spectrum, and lack of signiflcant variabil ity, the pulsar win d nebula of the 
Crab supernova remnant has been proposed as an X-ray standard fiux calibrator (jKirsch et al.l 120051 ). Observations 
over a large range of energies and with many diverse instruments reveal an intrinsic spectrum nearly consistent with 
a single power law; howeve r, the photon index and normalization determined by each detector exhibit small but 
not insignificant differences ()Kirsch et al.ll200"5l : iWeisskopf et al.ll2010f) . XMM-Newton EPIC-pn measurements, which 
require that observations are made in burst mode due to XMM-Newton's large collecting area, are best fit with a 
steeper tha n average photon index (F — 2.13 versus (F) = 2.08 for simultaneous fits to many instrument observations, 
iKirsch et~a l. 2005) that is driven by the shape of the spectrum from 0.5-2 keV. However, the residuals to this fit in 
their Figure 7 suggest that for energies above 2 keV, the XMM-Newton photon index is more in line with the average, 
and since we only consider energies above 2.3 keV, we need to determine what the Crab spectrum is measured to be 
in this range. 

Instead of fitting the Crab spectrum directly, we choose to compare the XMM-Newton data to Suzaku XISO data, 
which has been well calibrated and consistently fit over its energy range using observations of the Crab. Spectra from 
each ins trument ar e extracted from identical spatial regions at the center of the Coma cluster (specifically. Region 
10 from lWik et al.| [2009). where both the gas temperature and surface brightness are roughly constant. Fitting each 
spectrum from 2.3-12 keV with a single temperature APEC model, we find the XMM-A''ewto7i-derived temperature to be 
slightly, though not insignificantly, lower than the Suzaku temperature: 8.32 keV versus 8.90 keV. The lower EPIC-pn 
temperature is consistent with an increasingly larger effective area at higher energies relative to Suzaku^s high energy 
effective area; positing that XMM-Newton^s calibration is correct, the relative Suzaku effective area at larger energies 
should be increased as a function of energy, which would lower the fiux and therefore the temperature. We model 
this effect as multiplicative power law component to the Suzaku APEC fit, which simultaneously accounts for both 
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TABLE 3 
Swift BAT Error Factors 





Energies 


Flux Calib.'' 


„ b 
O^bgd 


PtSrc 


Thermal 


R25 


R30 


R35 


R40 


R45 


R60 


KW 


Band 


(keV) 


(cm~^) 


/m 




fm 








fm 








El 


14-20 


16.55 


2.37 


1.44 


2.09 


1.80 


1.77 


1.60 


1.73 


1.98 


2.08 


2.06 


2.22 


E2 


20-24 


11.28 


1.05 


1.32 


2.02 


1.75 


1.39 


1.53 


1.56 


1.22 


1.60 


1.57 


1.94 


E3 


24-35 


10.89 


1.63 


1.35 


1.97 


1.59 


1.41 


1.70 


1.82 


1.51 


1.81 


1.91 


2.08 


E4 


35-50 


10.46 


1.12 


1.31 


1.90 


1.66 


1.40 


1.59 


1.53 


1.18 


1.62 


1.64 


2.15 


E5 


50-75 


9.75 


1.02 


1.30 


1.88 


1.27 


1.59 


1.35 


1.54 


1.28 


1.45 


1.41 


2.04 


E6 


75-100 


13.00 


0.99 


1.29 


2.03 


1.34 


1.37 


1.52 


1.25 


1.36 


1.51 


1.57 


1.85 


E7 


100-150 


24.48 


1.75 


1.31 


2.67 


1.40 


1.47 


1.59 


1.57 


1.31 


1.52 


1.48 


1.72 


E8 


150-195 


75.19 


3.84 


1.31 


1.24 


1.42 


1.45 


1.26 


1.29 


1.23 


1.65 


1.56 


1.93 




Dittuse A'psE 


' (indcp. ot E): 




1.00 






4.56 


6.60 


8.98 


11.81 


14.94 


26.42 


41.03 



^ "Flux Calibration," defined sucfi that the incidence photon flux at the Earth (photons cm ^ s ^) is given by the 
BAT source count rate multiplied by Flux Calib. 
^ Units are 10~^ photons s~^ cm~^ 

the gradient and the overall cross-normalization between the two instruments. Fixing the APEC model parameters 
to those found with the fit to the XMM-Newton spectrum, we find this modification to the Suzaku calibration: 
f{Aeff) = 0.923(i?/l keV)°°*^. In other words, Suzaku spectra are fiatter than XMM-Newton spectra and have 
similar hard band fluxes; however, n ote that the XMM-Newton eff'ective area had previously been reduced by 15%, as 
per the analysis in iWik et all (|2OO90 . in order to match the 2-10 keV EPIC-pn and XISO fluxes. Dividing the XISO 
best-fit Crab spectrum of Ishida et al. [Fcrab,xiso(^^) = 9.51{E/lkeV)-'^ °^ , given in Suzaku Memo 2007-llE3| by 
f{Aeff) finally yields the correct parameters for the XMM-Newton fit to the Crab spectrum in the energy range of 
interest {E > 2 keV): 

/ ^ \ -2.095 

i^crab(£^) = 10.30 f^^j photons/cmVs. (Al) 

We take this equation as our canonical Crab spectrum. Then, the true flux of a source in each of the 8 bands is given 
by the BAT source count rate in that band divided by the obser ved Crab BAT coun t rate in that band, and multiplied 
by the spectrum in Equation (jAll) integrated over the band fsee lTueller et al.ll2010l eqns. 2-4). This conversion factor 
is reported in Table [Has the "Flux Calib.". Note that this photon index differs significantly from that assumed in 
iTueller et al.l (j2010( ). F = 2.15. Our goal here is to choose a photon index that agrees with the XMM-Newton data 
and best matches th e Cra b spectrum at the lowest BAT energies, where it is flatter than at harder energies (see the 
discussion in Section |A?2|) . 

While this method is the standard way to create spectra from BAT survey data, it is technically only valid for 
source spectra that have a shape similar to the Crab. Unfortunately, the BAT survey redistribution matrix, which 
could properly account for arbitrary spectral shapes, is more uncertain than the observed Crab fluxes (jTueller et al.l 
[20T0l). Thermal emission above ^10 keV is typically much steeper than the spectral index in Equation (IA1[) . so we 
would prefer to include an approximate redistribution matrix that will handle other such spectral models correctly. 
To do this, we take a standard response matrix for an on-axis source from a single observation and multiply the input 
energies by a smooth function so that the flux-converted Crab spectrum matches Equation (|A1|) to < 1% in all energy 
bands. The addition of this redistribution matrix has a minor effect on spectral fits generally, but it does improve the 
quality of fits using a thermal model, so we employ it throughout. 

A. 2. The Hard X-ray Spectrum of the Crab 

With this approach, accurate conversions from BAT count rates to true fluxes are not guaranteed. The goal instead 
is to match the BAT calibration with the EPIC-pn calibration, which ensures that spectral models can be applied 
seamlessly between the XMM-Newton and Swift spectra in joint fits. While fluxes quoted hereafter may differ from 
their true fluxes, the relative amounts of thermal versus nonthermal emission in the joint spectra - considering both 
their cross-normalization factor and shape - are carefully conserved. Ultimately, because our BAT calibration method 
relies on using the Crab as a flux standard, and since the true Crab spectrum is not known, the choice of a canonical 
Crab spectrum is at some level arbitrary. 

Even so, the hard band fluxes derived herein should be consistent wi th fluxes derived from other missions. Us- 
ing the same power law form as for Equation (jAll) . Kirs ch et al.l ()2005[) found a range of normalizations and pho- 
ton indices for several instruments that overlap with the BAT energy bands: Beppo-SAX PDS: 8.84£'~^-^^^; RXTE 
PCA: 11.02i;-2 i20. 7?XT£;HEXTE: 9.9£;-2 °90; INTEGRALISGKI: IbAlE-^-^^^; and INTEGRAL SPI: Ib.QE-^-^^^. 
These photon indices agree well wit h the shape of the total Crab spectrum determined by an earlier, detailed balloon- 
borne study, which found F = 2.12 (jBartlettlllQOl . Also, the Suzaku PIN fit of 10.93£'-2 090 jg consistent with both the 
scatter in the above results and our adopted spectrum. Ignoring the photon indices derived from XMM-Newton data, 
there seems to be a steepening in the Crab spectrum at higher energies, which means that our relatively fiat photon 
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index may over-predict harder band fluxes and thus enhance a potential nonthermal signal. Our Crab spectrum also 
has a slightly higher overall flux (in the 20-80 keV band), so that fluxes and upper limits may be biased high, though 
any such biasing would be well within the absolute calibration uncertainties of all the above missions. 

Assuming the INTEGRAL SPI observations are the most accura t e, due to its extensive ground calibration, we can 
compare our adopted model with the results of iJourdain fc Roaued ()2009f ). which show a clear steepening of the Crab 
spectrum above 100 keV. Below 100 keV, they flnd a photon index of 2.07, consistent with our value of 2.095. Their 
best-fit photon index above 100 keV is 2.24, which would lead our flux calibration to overestimate fluxes in our two 
highest energy bands. However, since we detect no flux in these bands, extended or other wise, even a large erroneous 
calibration - which is not the case here - would not bias or otherwise affect our results. IJourdain &: Rogues! ()2009f ) 
also fit a smoothly varying power law to their Crab spectrum. Not surprisingly, this model is sufficiently similar to our 
adopted Crab spectrum, particularly below 50 keV where the flux calibration is more important, that using it instead 
of our simple power law model would not lead to quantitatively different results. 



Very extended, diffuse emission is difficult to detect with coded mask instruments, since the shadow pattern of 
the mask on the detectors gets smeared out and the signal becomes indiscernible from the background. However, 
small scale extended emission is detectable, as long as its size is less than the minimum scale necessary to dilute the 
distinguishability of the mask pattern. In the following, we show that, by simulating extended sources as collections of 
point sources, this minimum scale is larger than our region of interest and that essentially 100% of the diffuse emission 
can be detected. 



For on-axis sources in the BAT FOV, the point spread function (PSF) has a full width half maximum (FWHM) of 
~22'. Because the sky image is basically the cross-correlation function of the coded mask with the count rates in the 
individual detectors, it does not represent the intensity per solid angle (i.e. within a pixel). Instead, a pixel value is 
proportional to the flux of a point source at that location. The width of the PSF is actually due to oversampling the 
sky plane, not the scattering of photons inside the instrument, and it depends on the size of individual mask element 
shadows on the detector relative to the size of detector pixels. 

As such, the PSF should not be summed in order to derive the source flux - this is provided by the central peak value 
- and its FWHM depends on the off-axis angle of the source. The distance between the mask and detector increases 
as the off-axis angle increases, so an angular separation at large off-axis angles produces a more dramatic shift in 
the shadow pattern across the detector pixels than more on-axis positions, which effectively reduces the oversampling 
factor and leads to narrower FWHM. Since survey images are created from many "random" individual pointings, each 
with a given source located at a different off-axis angle, the survey PSF will have an average FWHM and uniform 
shape, which is roughly Gaussian. Simple Gaussian fits to all the > lOcr sources in all 8 bands yield an average FWHM 
of 19f47, irrespectiv e of S/N or energy band, essentially identical to the value determined in the 22-month survey 
(|Tueller et al.ll2010l) . 

The PSF shape is described by a Gaussian to first order, which is not surprising given the non-repeating, randomly 
filled mask and the many pointings that contribute to the flux at each position. However, as is clear from the residuals 
to a Gaussian fit to the 14-20 keV band Crab data and another source (Cyg X-2) in Figure 13 deviations on the order 
of 1% of the flux exist (and are significantly larger than the root-mean square (RMS) of the background in this case). 
While this deviation does not strongly impact the flux of point sources, since only the maximum, central value maps 
to the flux, a diffuse source is composed of overlapping PSFs, where differences in the wings could affect the overall 
flux. The residual structure in the wings of the PSF is mainly eliminated by the addition of the two-part function: 



where r is the distance from the center in arcminutes, po and pi are fit parameters (the normalization and background, 
respectively), x = 2r/(1.19 FWHM) + 1 {x has units of radians), and ctpsf = FWHM/ (2^2 In 2). The improved fit for 
the Crab is illustrated in the very bottom left panel of Figure [Sj while in this case the fit is still not perfect, for other 
sources the fit is typically better (bottom right panel of Figure [3]). We take po to be our measurement of the flux. 
While the maximum of the PSF may not exactly correspond to po, since all fluxes are determined this way and are 
also related to the Crab fluxes, any such bias will cancel out during the conversion from BAT count rates to fluxes. 
Note that the most the additional terms to /(r) could affect a flux, assuming they, for some reason, poorly represented 
the true PSF shape, is at the < 1% level. 



As discussed in Section|31 we extract fluxes f or extended sou rces by fitting a priori model distributions, as opposed to 
using a method like the "CLEAN" algorithm (|Hogbomlll974[ ). which reconstructs fluxes from an unknown underlining 
distribution assuming the PSF shape only. The "CLEAN" method requires some flne-tuning, such as the region of 
extraction (for clearly detected sources, expanding the region-of-interest even a little beyond the wings of the source 
can significantly bias the derived fiux), "loop gain," and completion threshold. In our case, since there are only a few 



B. EXTRACTING BAT FLUXES FROM EXTENDED SOURCES 



B.l. BAT Point Spread Function 
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B.2. Tests of the Detection of Extended Sources 
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Fig. 9. — Fits to the BAT PSF. In the top and middle panels, the profiles of two point sources, the Crab and Cyg X-2, arc fit to a 
Gaussian. Large, regular residual s (m iddle panels) remain, which are mostly removed (at the < 1% level, bottom panels) after modifying 
our expression for the PSF (Eon. [BH . Lingering residuals, which are particularly significant for the Crab but only slightly noticeable for 
the other source, primarily result from spatial asymmetries due to the actual angular extent of the source, as in the case of the Crab, and 
the rectangular shape of the BAT instrument, which will cause oflf-axis sources to be more "squished" in one direction than another. In 
the latter case, the effect of summing many individual pointings with the detector in various orientations almost, but not entirely, removes 
this azimuthal component of the PSF shape. 



likely spatial distributions for the thermal and any potential nonthermal emission, we are less likely to produce biased 
fluxes by first assuming a spatial distribution than by using a method like " CLE AN." We represent a diffuse source 
as a collection of point sources, each of which is convolved by the PSF (Eqn. IB1[) and summed together. 

We now test whether diffuse sources are detectable over our scales of interest can be evaluated. In gen- 
eral, we treat extended emission as a collection of closely-spaced point sources, since existing software is built 
with these sources in mind. Point sources at any position in the BAT FOV are straightforward to simulate 
with the HEASOFT Swift task batmaskwtimg with the following options set: coord_type=tanxy; distance=le7; 
corrections=forward, unbalanced, flatfield; and rebalance=no. This task outputs the fraction of each detector 
pixel which is illuminated by the source given its position relative to the detector axis; a value of 0.45 means that 
55% of the detector area is shadowed by the mask. At this stage, the detector image can be multiplied by the counts 
or count rate of the source, and several such detector images representing different sources in the FOV can be added 
together along with a background - all including Poisson statistics. The background can then be fit and subtracted 
with the task batclean, and finally a reconstructed sky image can be produced via the task batf f timage. For now, 
to isolate the detectability of diffuse emission by the BAT, we simply add uniformly bright, perfectly known point 
source masks without background or source Poisson noise, to create circular, extended disks of various radii R. Images 
of the sky are constructed with batf f timage for each disk detector image, and the "observed" disk surface brightness 
profile is fit for as a function of radius. While even large disks {R > 10°) are visibly noticeable, sky images produced 
in this way include significant systematic effects induced by the large spatial extent of the emission, leading to high 
RMS noise that eventually destroys the disks' detectability. The recovered surface brightness, relative to the input 
level, of simulated diffuse disks of radius R are presented in Figure [TUl Error bars represent the simple error of the 
mean (ctrms/v'^; where N is the number of pixels used to determine i7rms)i and the smooth as opposed to random 
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Fig. 10. — Recovered surface brightness for simulated emission from a uniform surface brightness disk of a given radius. Error bars 
indicate the statistical error of the mean on the disk surface brightness, though the spatial fluctuations in the sky reconstruction behind 
this error are due entirely to systematic effects; the "noise" at any position in these simulations is determined by the fiux of all the other 
sources within the FOV, or in this case the other parts of the disk. The variation with radius is smooth instead of random due to the 
systematic origin of the fluctuations. For angular sizes of interest here, < 90', the intrinsic uncertainty in the recovered flux of an extended 
source, due to coded mask imaging techniques, is at most a few percent (based on the size of the error bars). Also, there is no loss in 
sensitivity to diffuse emission; all of the input flux is recovered, albeit with less and less precision for larger sources. 



variation around tlic input surface brightness results from their systematic nature. For sources in the size range of 
interest to us, R < 1?5, the intrinsic uncertainty due to the telescope design is < 3%. As the disk radius increases, the 
reconstructed surface brightness becomes less and less robust as there are effectively more sources (other parts of the 
disk) contributing systematic noise to a given location. Note, however, that all of the input source fiux is recovered. 

C. UNCERTAINTIES IN BAT FLUXES 

C.l. Flux Uncertainties for Point Sources 

The uncertainty in a given flux measurement is encoded in the RMS fluctuations in the local background 
(|Tueller et al.ll2010l ). These fluctuations represent both the statistical fluctuations from shot noise (dominated by 
the high background rate) and systematic error contributions from the sky reconstruction process. Due to the large 
number of individual pointings at nearly random positions, most systematic effects nearly average out and lead to a 
symmetric, nearly Gaussian distribution for blank sky regions. We calculate the RMS of the background (cTbgd) around 
Coma in an annulus of radius 15 < r < 100 pixels (42' < r < 4?67), as is typically done for sources in the BAT survey. 
The values of (Xbgd for each band are given in Table |31 While this annulus partially includes the region within which 
we are searching for a diffuse nonthermal signal, the lack of any obvious emission indicates that the derived errors 
could not be significantly biased. To ensure Cbgd is not biased by low level extended flux, we recalculated it inside an 
annulus of equal area with an inner radius of 90' and found a nearly identical value of Ubgd in all 8 bands. 
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C.2. Flux Uncertainties for Extended Sources 

The error for a point source, or more correctly the error in the value of a given pixel, presented in Section IC.ll 
does not directly apply to extended sources. Also, we cannot take the standard error from spatial fits, using the 
point source error as the error for the flux in each pixel, because nearby pixels are correlated. Helpfully, the expected 
error for diffuse sources has already been derived bv lRenaud et al.l (|2006f) for the IBIS coded mask instrument onboard 
INTEGRAL. In their appendices, they derive source fluxes and errors in reconstructed sky images from detector images 
and find the straightforward result that the error in a flux measurement from a n extended source is proportional to 
its spatial area normalized by the area of the PSF function (jRenaud et al.|[2006l Eqn. B3). Specifically, 



Ccxt — fbffd 



V^Vps^, (CI) 

where iVpgp is the area of the source normalized by the PSF area. If Ix{^) is the surface brightness of the source 
convolved with the PSF, I^^^ is its maximum value, ft is the solid angle, and f{^) is the PSF, then 

_ Jixin)dn 

Before we generally apply Equation (|C1|) to Swift BAT data, we test whether this prescription does in fact apply 
to extended sources in the Swift BAT. For each of the diffuse thermal (Section 13.1.11) or nonthermal (Section I3.1.2P 
spatial models described above, we created 121 simulations of the extended source and a number of other point sources 
in the FOV, each with detector counts for the source fluxes and background randomly assigned (taken from a Poisson 
distribution). Each simulated observation is made from a unique position on an 11 x 11 grid, with the relative positions 
of all the sources kept intact. Both the variation of position relative to the telescope axis and the inclusion of point 
sources are necessary to fully recreate the systematic contribution to the error. The simulated detector images are then 
background subtracted and converted into sky images via the procedure outlined in Section [B.2l The total flux of the 
diffuse sources is chosen so that the signal-to-noise ratio is ~ 20 — 50. To check Equation (|Cip . fluxes of all the sources 
are measured in each simulated sky image, and the average standard deviation of the point source fluxes are compared 
to the standard deviation of the diffuse source flux. We find that the estimated errors for the diffuse models (Thermal 
band El, Disks R25-R60, and KW) generally fall below the expected trend with iVpsF in Figure [TTJ This discrepancy 
may be due to the number of simulations we were computationally limited to per forming - the distribu tion of fluxes 
is only roughly Gaussian - or it may represent a true deviation from the results of iRenaud et al.l ()2006D . However, to 
be safe we use Equation (jC3|) to calculate the error of fluxes extracted with the corresponding model. 

There is one additional modiflcation to errors on fluxes extracted with our methodology. Because we flt a spatial 
model, convolved by the PSF, to the BAT image data, the error in the flux is not just the standard deviation of 
nearby background pixels, but it depends on how the model is flt to all the pixels. For example, the distribution of 
normalizations from many fits of a Gaussian function to random data of mean zero and standard deviation Cbgd will 
not equal CTbgd but some value < <7bgd depending on the pixel scale. A delta function, or Gaussian of width zero, will 
produce a distribution consistent with tTbgd, since this is identical to measuring the standard deviation, but anything 
wider finds an average over several pixels, and therefore the distribution of normalizations will tend to be closer to the 
mean of the random pixels. For our purpose, where the normalization is related to the source flux, the correct error of 
a fiux should come from the distribution of model fits to background (empty) regions of the survey, which may not be 
equivalent to Ubgd- Unlike in the above example, neighboring pixels in the survey are correlated due to oversampling 
- this is essentially the origin of the PSF - and so the standard deviation of model normalizations will be affected by 
this correlation. Generally, the distribution of normalizations will be larger than Ubgd in this case, as minimization 
will be more influenced by the larger fluctuations near the flux extraction region. The net effect does not significantly 
change the error distribution shape, but simply inflates the effective standard deviation by some factor, /„, which is 
both model-dependent (varying from 1.4tTbgd for a point source to 2.24(Tbgd for the KW model) and energy dependent 
since the noise properties vary slig htly from band-to-band. The total flux uncertainty for a diffuse source in the BAT 
survey is adjusted from Equation (|C1I) to become 



^diffuse — /mfcxt — /mCbgd V^VpsF . (C3) 

The precise value of /,„ is determined from the standard deviation of flts to 100 blank sky regions, in which we avoid 
obvious (> 5a) sources and the Galactic plane {b > 20°). These factors are reported for each band in Tabled Not 
including this factor in the flux error estimate results in spectral flts with unacceptably high values. 
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